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Abstract 

We analyze the possibility to measure small thermal effects in the Casimir force between metal 
test bodies in configurations of a sphere above a plate and two parallel plates. For sphere-plate 
geometry used in many experiments we investigate the applicability of the proximity force approx- 
imation (PFA) to calculate thermal effects in the Casimir force and its gradient. It is shown that 
for real metals the two formulations of the PFA used in the literature lead to relative differences 
in the obtained results being less than a small parameter equal to the ratio of separation distance 
to sphere radius. For ideal metals the PFA results for the thermal correction are obtained and 
compared with available exact results. It is emphasized that in the experimental region in the 
zeroth order of the small parameter mentioned above the thermal Casimir force and its gradient 
calculated using the PFA (and thermal corrections in their own right) coincide with respective 
exact results. For real metals available exact results are outside the application region of the PFA. 
However, the exact results are shown to converge to the PFA results when the small parameter goes 
down to the experimental values. We arrive at the conclusion that large thermal effects predicted 
by the Drude model approach, if existing at all, could be measured in both static and dynamic ex- 
periments in sphere-plate and plate-plate configurations. As to the small thermal effects predicted 
by the plasma model approach, the static experiment in the configuration of two parallel plates is 
found to be the best for its observation. 

PACS numbers: 31.30.jh, 12.20.Ds, 12.20.Fv, 42.50.Nn 
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I. INTRODUCTION 

The Casimir effect [l| is now universally known as one of the most extensively studied 
manifestations of zero-point oscillations. It attracts considerable attention in quantum field 
theory, gravitation and cosmology, atomic physics and optics, condensed matter physics, 
and in applications to nanotechnology (see monographs (2|-|6|). The fundamental theory 
describing both the van der Waals and Casimir forces between two semispaces with planar 
boundaries was developed by Lifshitz j?, 8|. At present the Lifshitz theory is generalized for 
material bodies with curved boundary surfaces (see, e.g., Refs. jo], Q). Over a long period 
of years only a very limited experimental information about the Casimir force was available. 
Recently, however, a significant advance has been made in the experimental study of this 
phenomenon reflected in review IjJ. In principle, the theoretical and experimental progress 



taken together allows computation and measurement of Casimir forces between bodies of 
complicated geometrical shape made of different materials. In order for such computation 
to be made, one should know the reflection amplitudes on the boundary surfaces at all 
Matsubara frequencies. At flrst sight this should present no problems because the reflection 
amplitudes can be either measured directly or calculated using measured material properties, 
such as frequency-dependent dielectric permittivities. At this point, however, difficulties 
emerge unexpectedly which are connected with the role of relaxation properties of charge 
carriers in the Casimir effect. 

Beginning in 2000, the influence of relaxation on the thermal Casimir force was hotly 

n 

debated. Bostrom and Sernelius [Iz] have noticed that the substitution of the dielectric 
permittivity of the Drude model with nonzero relaxation parameter into the Lifshitz formula 
results in a decreasing magnitude of the Casimir free energy and force as a function of 
temperature over a wide region of separations. This is in contradiction with the case of ideal 
metal plates or plate materials described by the nondissipative plasma model, where the 
magnitudes of the Casimir free energy and force are monotonously increasing functions of the 
temperature 3, Q]- At large separations (or in high temperature limit) the magnitudes of 
the Casimir free energy and force between metal plates calculated using the Drude model are 
by a factor of 2 less then the same quantities calculated for ideal metals or metals described by 
the plasma model. It was also shown that for metals with perfect crystal lattices the Lifshitz 
theory combined with the Drude model violates the Nernst heat theorem [g], [ll|, [is], [l^ . For 
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metals with impurities leading to a nonzero relaxation at zero temperature Nernst's theorem 



was shown to be followed 
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18| . This, however, does not solve the problem because perfect 



crystal metal is an idealized model of a truly equilibrium system with nondegenerate ground 
state for which the laws of thermodynamics must be satisfied. 

A large thermal effect in the Casimir force between metals at separations of a few hundred 
nanometers predicted by the Lifshitz theory combined with the Drude model was experi- 



mentally excluded 



19| by indirect dynamic measurements of the Casimir pressure in the 



configuration of a sphere above a plate of a micromachined oscillator. Later this experiment 
was repeated for two more times with increased precision. The exclusion of the thermal 
effect predicted by the Drude model was confirmed at a 95% confidence level 2o|, 21] and at 



a 99.9% confidence level [22|, |23|. The same measurement data were found to be consistent 



with the Lifshitz theory combined with the plasma model 
that the comparison of experiment with theory in Refs. 
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t is important to keep in mind 



23| was based on the use of an 



approximate method, the so-called proximity force approximation (PEA) 



ay 



24| because 



at that time for sphere-plate configuration an exact theory was not available. At the moment 



there is a conceptual possibility to compute the thermal 
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Casimir force between a sphere and 



261], but the region of experimental 



a plate made of real metals with no use of the PFA 
parameters is not yet achieved due to computational difficulties. 

When it is considered that the Drude model correctly describes the relaxation of con- 
duction electrons at low frequencies, the contradiction with basic laws of thermodynamics 
and disagreement with the experimental data outlined above are puzzling. These problems 
were dramatized by the demonstration that the inclusion of dc conductivity of dielectric (or 



dielectric-type semiconductor) materials into the model of die 



theory also results in a violation of the Nernst theorem 



ectric response in the Lifshitz 



27H30|. From the experimental 



side, it was shown that the measurement data of the experiment on optical modulation of 



the Casimir force between Au sphere and Si plate with light |31l . 1321 ] excludes the Lifshitz 
theory taking into account the dc conductivity of dielectrics at a 95% confidence level. A 
similar result was obtained from the measurement of the Casimir-Polder force between the 
Bose-Einstein condensate of ®^Rb atoms and Si02 plate 33|. Here, the Lifshitz theory tak- 
ing dc conductivity of Si02 into account was experimentally excluded at a 70% confidence 
level ^]. It is pertinent to note that while the comparison of the optical modulation ex- 
periment with theory uses the PFA, the measurement results for the Casimir-Polder force 
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were compared with the exact Lifshitz formula for atom-wall interaction. One can summa- 
rize that experiments with metals, semiconductors and dielectrics exclude the influence of 
dissipation of conduction electrons on the Casimir force (the statement on the opposite in 



the Introduction to Ref. 



26| is a typo). 



The conflict between such a fundamental theory, as the Lifshitz theory, thermodynamics 
and the experimental data of several experiments is a problem of great concern. Because of 
this, a lot of attempts to resolve this problem has been undertaken. Specifically, it was even 
suggested 35|, |36| to modify the Lifshitz theory by including into consideration screening 



effects and diffusion currents. It was noted, however, that the modifications proposed do not 
alleviate contradictions with thermodynamics and the experimental data (discussion on this 



subject can be found in Refs. 



37H42|). On the other hand, it was suggested j43| to modify the 



Planck distribution law by taking into account "saturation effects" . There were also attempts 
to soften contradictions with thermodynamics by reformulating the problem 44j | an d by 



finding additional statistical arguments in favor of the Drude model 45 



46|. In Ref. 



in addition to the usually used exponential screening, the so-called algebraic screening in 
atom- wall interaction was considered. As a result, linear in temperature thermal correction 
to the Casimir-Polder force at short separations was predicted similar to that predicted by 
the Drude model. Note that the algebraic screening is connected with nonanalytic terms 
in the small wavenumber expansion of the dielectric permittivity. It was finally suggested 
48l | that for the resolution of the problem some concepts of statistical physics related to 
the theoretical description of the interaction of classical and quantum fluctuating fields with 
matter might need a reconsideration. 

The prospects for a pure theoretical resolution of the above problems seem dim at the 
moment. In this situation any additional experimental evidence could be very useful. In this 
paper we analyze the possibility to measure thermal effects in the Casimir force on the basis 
of already created and used experimental setups. We stress that in the experiments 0- 
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31 



321] mentioned above the large thermal correction to the Casimir force, as predicted by 



the Drude model, was excluded. However, these experiments were not of sufficient precision 
to measure the thermal effect for metals predicted by the plasma model or for dielectrics with 
dc conductivity omitted (till the moment the thermal effect was measured in the Casimir- 
Polder force alone [3^]). Keeping in mind that there is some confusion in the literature 
concerning the use of the PFA, we present two (not equivalent) formulations of the PFA and 
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clarify which of them was really used in the comparison between experiment and theory. 
We especially analyze the calculation results for the thermal correction to the Casimir force 
in sphere-plate configuration found using the PEA and explain when they are meaningful. 
The obtained conclusions are confirmed by the comparison with exact results for the thermal 
contribution to the Casimir force between a sphere and a plate. We show that at the moment 
it is not possible to measure small thermal effects in the Casimir force (or its gradient) in 
sphere-plate geometry, as predicted by the Lifshitz formula combined with the plasma model. 
We also discuss the configuration of two parallel plates in both dynamic and static regimes. 
According to our results, small thermal corrections to the gradient of the Casimir pressure 
in the dynamic regime is suppressed. The only way to measure small thermal effects in the 
Casimir pressure is suggested by the configuration of two parallel plates in the static regime. 

The paper is organized as follows. In Sec. II two different formulations of the PFA 
are considered and applied to the Casimir force in sphere-plate configuration. Section III 
discusses the same subject with respect to the gradient of the Casimir force between a sphere 
and a plate. In Sec. IV the relationship between the PFA and the exact results is presented 
for a sphere and a plate made of ideal metal. The applicability of the PFA to describe 
thermal corrections to the Casimir force between a sphere and a plate made of real metals 
and the possibility to measure the thermal effect are considered in Sec. V. In Sec. VI we show 
that the static Casimir configuration of two parallel plates is preferential for the observation 
of a small thermal effect in the Casimir pressure. Section VII contains our conclusions and 
discussion. 



II. THE PROXIMITY FORCE APPROXIMATION FOR THE CASIMIR FORCE 
BETWEEN A SPHERE AND A PLATE 

Different authors vary somewhat in the meaning of the term "PFA". In fact the term 



"proximity force theorem" (later changed for PFA) was introduced in Ref. [24| where the 



so-called Derjaguin method [49| was applied in order to calculate the force acting between 
curved surfaces by using the known force per unit area of plane parallel plates. In the Der- 
jaguin method, the unknown force between the elements of curved surfaces is approximately 
replaced with a known force per unit area of plane surfaces at respective separations. In 
application to a sphere of radius R above a plane surface of a plate z = the Derjaguin 
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method represents the force between them in the form 

F,p(a,T) = [ daP{z,T). (1) 



Here, da is the element of plate area, S is the projection of the sphere onto the plate, a is the 
shortest separation between the sphere and the plate, and P{z,T) is the force per unit area 
of two plane parallel plates at a separation z at temperature T (i.e., the Casimir pressure). 
Choosing the origin of a cylindrical coordinate system on the plane z = under the sphere 
center, the coordinate z of any point on the sphere is given hy z = R + a — {R^ — p^Y^"^- 
Then Eq. ([1]) leads to 

F,p{a,T) = 271 / pdpP{z,T) (2) 
Jo 



pR+a 

2n / {R + a- z)P{z,T)dz. 

J a 



Keeping in mind that the thermal Casimir pressure is connected with the free energy per 
unit area of two parallel plates as 

P{z,T) = -^^^^jf^, (3) 
and integrating by parts in Eq. ([2]), one arrives at 

pR+a 

Fsp(a,r) = 27ri?J-pp(a,T) - 2n dzJ^^^{z,T). (4) 

J a 



This generalizes Eq. (20) of Ref. 50|] related to the nonretarded case. 

Further simplification of Eq. (jlj) can be achieved when it is assumed that the free energy 
J^pp(2;, T) is a quickly decreasing function of z and drops to zero on the characteristic length 
of about the sphere radius R. Let us consider the case of the free energy of the Casimir 
interaction given by the Lifshitz formula, 

:f,,{z,t) = / ^aEi^(i-^'^'"''^)- (5) 

Here, fee is the Boltzmann constant, the prime near the summation sign multiplies the 
term with / = by 1/2, fc_L is the projection of the wave vector on the plane of plates, 
and qi = {kl + ^f/c^Y^^ where = 2nkBTl/h with / = 0, 1, 2, . . . are the Matsubara 
frequencies. The reflection coefficients for the two polarizations of the electromagnetic 



6 



field, transverse magnetic {a = TM) and transverse electric {a = TE), are expressed in 
terms of the dielectric permittivity along the imaginary frequencies, 



qi - h 



(6) 



k, 



kl + e{i^i) 



if 



2^ 1/2 



In order to simplify Eq. (|4]), we use an expansion in power series in Eq. 



oo oo 



J^pp(^,T) = 
Then from Eq. ([7]) one finds 



2tt ^ — ' — ' n 

1=0 n=l 



f-R+a h T ^ I ^ 1 



(7) 



oo oo 



X 



k^dkj 

qi 



1=0 n=l 



(8) 



When R goes to infinity with a and T fixed, the contribution to /(a, T) of the second term 
in square brackets on the right-hand side of Eq. ([8]) vanishes as {R + a)~^. This means that 
for large R the value of I{a,T) is determined by the first term in square brackets and can 
be considered as independent on R. Hence, rewriting Eq. in the form 

F,p(a, T) = 27ri?J-pp(a, T) [l + /^""f ^, ] , (9) 

one concludes that in the limit of large R it holds I{a,T)/ RJ^pp{a,T) C / R where C is 
some constant. 

Let us now consider the behavior of the quantity I/RJ^pp in the limiting case a — )■ 
keeping R fixed. This is a nonrelativistic limit where 

m 

H 



with the Hamaker constant defined by 



I{a,T) 



RH 

Uira^R + a) 



(10) 



I dil y^dy 



-1 



(11) 
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From Eq. (fTO|) it follows I{a,T)/RJ^pp{a,T) ~ Ca when a vanishes. Thus Eq. (|9]) can be 
rewritten as 



Fsp(a,r) = 27r/?J-pp(a,T) ^l + /(a,T)-J , (12) 

where /(a, T) is scarcely affected by the sphere radius R. For all models of dielectric permit- 
tivity used to describe metals the magnitude of f{a,T) is less than unity over wide ranges 
of experimental parameters (see below the results of numerical computations). 

Under the condition a <^ R which is usually valid in experiments on measuring the 
Casimir force one may neglect the term of order a/R on the right-hand side of Eq. f|T2l) and 
arrives at 

Fsp(a,T) =27ri?J-pp(a,T). (13) 



Just this equation was called the "proximity force theorem" in Ref. 2^ and heavily used in 



the comparison of experiment with t 
plate geometry (see, e.g., Refs. [31 



leory m a. 



32 



52 



1 measurements of the Casimir force in sphere- 



62| ) . As can be seen from the above derivation. 



Eq. f|T3l) follows from Eq. ([T]) under some conditions, but is not equivalent to it. Because 
of this, to avoid confusion, Ref. 63| suggested to call (II]) the most general formulation of 
the PFA and ( |T3l) the simplified formulation of the PFA. Keeping in mind that the used in 
experiments simplified formulation is obtained by disregarding contributions of order a/R, 
it would be meaningless to attribute physical meaning to any terms of order a/R in the 
thermal Casimir force Fgp calculated by using Eq. f|T3l) (see discussion in Sees. IV and V). 



64 



65| is connected with the 



One further version of the PFA discussed in the literature 
choice of parallel surface elements in Eq. ([T]). In the original Derjaguin method [49| used 
by us the surface elements representing the sphere are parallel to the surface of the plate 
z = 0. In this case S is a part of the plane z = (the so-called plate-based PFA). If, 



however, the lower half of the sphere is chosen as S [6J], the plane elements representing the 
sphere are tangential to it and respective twin elements of the plate are tilted by different 
angles with respect to the plane z = (the so-called sphere-based PFA). In both cases, the 
distance between the elements is measured along the normal to S. It was shown 64] that 
both the plate-based and sphere-based PFA lead to coinciding results in the zeroth order 
of a/R, i.e., to Eq. f|T3l) . The form of the function f(a, T) in Eq. f|T2|) for both versions of 
the PFA is, however, different. Basing on this, Refs. 65|, |66|| considered the results for the 
Casimir energies and forces obtained by PFA as ambiguous. The differences between the two 
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versions of the PFA were treated as some "error bars" inherent to this approximate method. 
Keeping in mind, however, that the comparison of experiment with theory in sphere-plate 
geometry is based not on the general formulation of the PFA ([T]), but on an unambiguous 
simplified formulation f[T^ . the discussion of inherent to PFA errors is immaterial. In fact, 
when speaking about the Casimir force, only the zeroth order m a/R results in any of the 
PFA formulations are of physical significance and only under the condition a <^ R. 

The situation changes drastically when, instead to the Casimir force, the PFA is applied, 



for instance, to the gravitational force. According to Refs. [63|, |67l], the most general for- 
mulation of the PFA ([T]) leads to an exact result for the force between a sphere and a plate 
for all conservative volumetric forces, particularly for the gravitational force. In so doing, 
however, only the original Derjaguin choice of S (plate-based) must be used. This makes 
the plate-based version of the PFA preferable in comparison with the sphere-based version. 

Now we present the results of numerical computations for the function /(a, T) defined in 
Eq. ( fT2|) in typical regions of experimental parameters. Computations were performed by 
using the Lifshitz formula ([5]) and Eq. ([8]) for a sphere and a plate made of Au. The dielectric 
properties of Au were described by using three different models. As a crude approximation, 
the model of ideal metal bodies was used leading to rTM(^6?^±) = 1; ''^te(^6?^±) = ~1- 
A frequently used description obtains the dielectric permittivity e{i^i) by means of the 
Kramers-Kronig relation 

.«0 = 1 + -/ -^J^d-^ (14) 

where lme{u) is taken from tables of the optical data for Au 681] extrapolated to low 
freuqencies by means of the Drude model 

2 

Imeico) = ""f (15) 
Here, the plasma frequency and the relaxation parameter of Au are given by = 9.0 eV, 



7 = 0.035 eV 69|. As one more alternative description used in the literature we have applied 
the generalized plasma-like model 

^m-^4^t':jr^,, (16) 

where lOj ^ are the resonant frequencies of the oscillators describing core electrons, 7j 
are the relaxation frequencies, and gj are the oscillator strengths. The values of all these 
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parameters for Au can be found in j^, [2^. Note that there are proposals in the hterature 
for alternative dielectric functions taking into account the effect of spatial nonlocality (see, 



for instance, Refs. 142 



70 



7l|). However, for metallic test bodies these dielectric functions 



41 



70 



7l|. Because 



predict precisely the same thermal effect as the Drude model approach 
of this, they do not require a special consideration here. 

In Fig. 1 the computational results for f{a, T) are shown as a function of (a) separation at 
room temperature T = 300 K and (b) temperature at a separation a = 2 fim. In both cases a 
sphere with R = 100 /im radius is used. The dotted, dashed and solid lines demonstrate the 
results obtained using ideal metals, Au described by the optical data extrapolated by the 
Drude model and Au described by the generalized plasma-like model, respectively. As can be 
seen in Fig. 1(a), in the region of separations from 100 nm to 5 /im it holds 0.5 < \ f{a,T)\ < 1. 
It can be easily verified that in the limiting case a — )■ the function f{a,T) — )■ —1/2 for 
ideal metals (the dotted line) and f{a,T) — )■ —1 for real Au independently of the model 
used for its description (the dashed and solid lines). From Fig. 1(b) it is seen that at any 
temperature from absolute zero to 300 K the magnitudes of /(a,T) remain less than unity 
and the difference between various models of dielectric properties disappears with vanishing 
temperature. 

Now we consider the measure of dependence of the function /(a, T) on the sphere radius. 
The computational results for /(a, T) as a function of logiQ{a/R) are shown in Fig. 2 as the 
three groups of lines (dotted, dashed and solid for ideal metals and real Au described using 
the Drude and plasma models as explained above) numerated 1, 2, and 3 for fixed separations 
a = 0.1, 2, and 5 fim, respectively. Computations are done at room temperature T = 300 K 
for a/R varying from 10~^ to 0.05. For example, as is seen in Fig. 2, at a = 0.1 /im (the 
group 1 of lines) there is almost no dependence of /(a,T) on R for radii larger than 2 /im. 
For a = 2 /im (the group 2 of lines) f{a,T) does not depend on R for radii R > 200 /xm. 



We remind that typical sphere radii are R = 100 /xm in the experiments 31 



32 



53 



62| and 



R = 150 /im in the experiments 20l-l23l|. performed at separations from less than 100 nm 
to a few hundred nanometers. In Fig. 2 it is seen also that the relative magnitude of the 
function /(a,T) computed using different models of dielectric permittivity of Au depends 
on separation. For example, at a = 0.1 /im \f{a,T)\ computed using the model of ideal 
metals (the dotted line) is less than using the generalized plasma-like model (the solid line) 
and using the tabulated optical data extrapolated by the Drude model (the dashed line). 
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At separations a = 2 and 5 /im the magnitude of /(a, T) computed using the model of ideal 
metals is sandwiched between lines computed using the two models of real Au. 

The size of possible corrections to the simplified formulation of the PEA (IT^ has been 
investigated experimentally by measuring the Casimir force between an Au-coated plate 
and five Au-coated spheres with different radii using a micromachined oscillator 72|. In so 
doing spheres with radii R = 10.5, 31.4, 52.3, 102.8 and 148.2 fim were used. The obtained 
constraint \f{a,T)\ < 0.6 for a < 300 nm, T = 300 K is in very good agreement with the 
computational results in Fig. 1(a). Thus, the experimental data are in favor of the simplified 
formulation of the PFA. 



III. THE PROXIMITY FORCE APPROXIMATION FOR THE GRADIENT OF 
THE CASIMIR FORCE BETWEEN A SPHERE AND A PLATE 

In many experiments using the sphere-plate configuration separation distance between 
the test bodies was varied harmonically (see, e.g., Refs. , Il9l-l23[|). In this case not 



the Casimir force but its gradient is the physical quantity immediately connected with the 
frequency shift of a micromachined oscillator. The gradient of the Casimir force acting 
between a sphere and a plate can be found by differentiating Eq. (jlj) with respect to a and 
taking into account Eq. ([3]), 



dF,p{a,T) 



da 

Using Eq. this can be rearranged to the form 

^^y^ = -2vri?P(a,T) 
da 

where the following notation is introduced: 



27rRP{a, T) + 27rJpp(a, T) - 27rJpp(i? + a, T). (17) 



J(a,T) 



RP{a,T) 



(18) 



, rp OO OO 1 poo 



r 

a 



1=0 n=l 

X |^e~^'^'"" — e"^'^'"'-'^"'""^] . (19) 

Similar to Sec. II, it can easily be shown that at large i?, with a and T fixed, it holds 
J{a,T)/ RP{a,T) C / R. On the other hand, in the nonrelativistic limit 

H ^, RH(R + 2a) , , 

11 



where H is defined in Eq. (fTTj) . From this it follows that J{a,T)/RP{a,T) ~ Ca when a 
vanishes and R is kept constant. Thus, Eq. f lTSj) can be rewritten in an equivalent form 



dF,p{a,T) 



a 



2nRPia,T) l+p(a,T)-J, (21) 



da 

where p{a, T) is scarcely affected by R. Below we demonstrate that in wide ranges of exper- 
imental parameters p{a,T) is a very slowly varying function and \p{a,T)\ < 1/2. Because 
of this, under the experimental condition a <^ i? we can neglect the term of order a/R in 
Eq. (pT]) and arrive at the equality 

^^'"q""^^^ = -27rRP{a, T), (22) 

which is the simplified formulation of the PFA for a gradient of the Casimir force in sphere- 
plate configuration. Equation (1221) was used for the comparison of experiment with theory 



in dynamic measurements of the Casimir pressure [19H231. 17 

We have performed numerical computations of the quantity p{a, T) as a function of sepa- 
tion and temperature for a sphere of fixed radius R = 100 fim. In Fig. 3(a) the computational 
results for p{a, T) as a function of a are presented at T = 300 K. Figure 3(b) shows p(a, T) as 
a function of T at a = 2 /xm. In both cases dotted, dashed and solid lines indicate the use of 
ideal metal surfaces and Au surfaces described by the Drude and plasma model approaches. 
As can be seen in Fig. 3(a), within the separation region from 100 nm to 5 /im the function 
|p(a, T)| varies between 0.31 and 0.48. When the temperature increases from absolute zero 
to T = 300 K, |p(a, T)| remains sandwiched between 0.31 and 0.39 [see Fig. 3(b)]. This 
demonstrates that possible corrections to the simplified formulation of the PFA (!22|) do not 
exceed a/R. 

It can be easily seen that the magnitude of p(a, T) almost does not depend on the radius 
of the sphere. In Fig. 4 we present the computational results for p(a, T) as a function of 
log|o(ct/-R) at T = 300 K. The groups of lines numbered 1 and 2 are for the separations 
a = 0.1 and 5 /im, respectively. The dotted, dashed and solid lines have the same meaning 
as in Figs. 1-3. As is seen in Fig. 4, for the group of lines 1 there is no dependence on a/R 
over the whole range from 10~^ to 0.05. Thus, at a = 0.1 /im, p(a, T) does not depend on R 
for R>2 /im. For the group of lines 2, p{a, T) does not depend on R for 10~^ "£ a/R < 10~^. 
At a = 5 /im this leads to R> 5000 /im. From Fig. 4 it also follows that relative magnitudes 
of the function p{a, T) computed using different models of dielectric properties of metal 
depend on separation. 
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The experimental constraint on the magnitude of the correction to Eq. (l22l) was obtained 



in Ref. 



72| by using several spheres with different radii (see Sec. II). In the separation 



region a < 300 nm at T = 300 K it was shown that \p{a,T)\ < 0.4 at a 95% confidence level. 
This is in very good agreement with the computational results in Fig. 4 and provides the 
experimental confirmation for the simplified formulation of the PFA in application to the 
gradient of the thermal Casimir force between a sphere and a plate. 



IV. RELATIONSHIP BETWEEN THE PROXIMITY FORCE APPROXIMATION 
AND EXACT RESULTS FOR IDEAL METALS 

In what follows we apply the PFA to calculate the thermal Casimir force between a sphere 
and a plate made of ideal metal and discuss the possibility to describe the thermal correction 
using this approximate method. Using the notation ([HD, Eq. (jlj) for the thermal Casimir 
force can be presented in the form 

Fsp(a, T) = 27ri?J-pp(a, T) + 27r/(a, T) (23) 
= 27r/?Jpp(a, T) + 27rX{a, T) - 2r(X{R + a,T). 

For ideal metals, using the dimensionless variables 

y = 2aqu ^ = rJ (24) 
c 

in Eq. ([8]), the quantity X(a, T) is given by 

X(a,T) = ^E'f:^ rdye--y 
47ra frf ^ Jr„i 



1=0 n=l 
oo 



I. 1^ °° 

l=—oo ' 



or the transition to 



where Lin{z) is a polylogarithm function. 

Equation fl25l) can be presented in an equivalent form convenient 
the low-temperature limit using the Poisson summation formula js], |6|. According to this 
formula, if c{a) is the Fourier transform of a function b{x), i.e. 

c{a) = ^ [ b{x) e-'""dx, (26) 
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then it follows that 

oo oo 

b{l) = 27rJ2 c(27r/). 

l = — OD l = — 00 

Now we return to Eq. (125 p and put 



m 



Jra\l\ 



Keeping in mind that h{l) = b{—l) and using Eq. fl26l) . one obtains 

coo 



c{27tI) 



he 



327r3a2 



b{x) cos(27r/x) dx 

POO 

dv cos{ltav) / dy Li2{e~^). 



(27) 



(28) 



(29) 



Here, we have introduced a new variable v = TaX and the notation ta = Ta/T, where 
ksTa = he/ {2a) is the effective temperature related to the separation distance between the 
sphere and the plate. Note that in terms of this notation it holds Tq = 211 /ta = 2TtT/Ta. 
Taking into account that the quantity c(27r/) is an even function of its argument and using 
Eqs. ( 1271) and ( l29|) . we arrive at 

oo 

X(a,T) = 47r^'c(27r/) (30) 

1=0 

he ' 

y I dy Li2{e~^) / dv cos(/taf), 



where we exchanged the order of integrations with respect to v and to y. 
Equation ( !30l) can be represented in the form 

he 



where 



X(a,T) 



^0 



\xo + Y.Xi 



1=1 



ydyU2{e ^), 



(31) 



(32) 







/"OO 

It. 



'^a JO 



dyU2{e sm{ltay). 



Direct calculation leads to 



n=l ''^ 



xdxe ^ = — , 
90' 



(33) 



1 °° 1 
Ita ^ 

n=l 



dxe ^ sin 



It^x 



n 



E 



TT 



^^,n\n^ + lHl) 21X QlHl 



7r coth(7rZt(i) 
2 Ptl ■ 
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where x = ny. Substituting Eq. ( 133|) into Eq. ( I3T]) . we obtain 



X{a,T) 



1440a2 



-I 



(34) 



90 coth(7r/ta 



ttHI 

1=1 



+ 



In the low-temperature limit it holds T <^ T^, ta ^ 1 and Eq. f lM|) results in 



X{a,T) 



TT^hc 



1440a2 
90C(3) 



1 + 5 



(35) 



where ({z) is the Riemann zeta function. 

The quantity X{R + a, T) in Eq. (12^ is given by Eq. 
ta is replaced with = Tr/T and Tr is defined from 

he 
2R 



T 



where a is replaced with R + a, 



ksTR 



(36) 



2{R + a) 

The sphere radius introduces a second characteristic temperature into the problem which is 
much less than in experimentally relevant situations. As an example, for typical spheres 
used in experiments on measuring the Casimir force R = 100 /im and a = 100 nm resulting 
in Tr = 11.4 K and Ta = 11400 K. Thus, for experiments performed at room temperature, 
T = 300 K, it holds T <^ Ta, but T ^ Tr, i.e., a high-temperature regime with respect to 
the sphere radius, but a low-temperature regime with respect to the separation between the 
sphere and the plate. 

Nevertheless, we begin with the case of extremely low temperature with respect to both 
parameters, T ^ Tr ^ Ta, which is achieved only well below 1 K. In this case the quantity 
X{R + a, T) is also given by Eq. ( I35l) where a is replaced with R + a and Ta is replaced with 
Tr. We also take into account that under the condition T <^ Ta the free energy between 
two plane parallel plates in Eq. (!23|) is given by js, 6| 



•?>p(a,7') 



45C(3) 



T 

Ya 



(37) 



Substituting Eq. (1371) and Eq. (1351) (with the characteristic temperatures Ta and Tr) into 
Eq. fl23|) . we arrive at the result 



i^sp(a,T) 



TT^hcR 
360a3 



a a 
2R ^ 



T 



(38) 
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Here, we have included the main T- independent term of order a/R and the first nonvanishing 
term depending on temperature. The latter is of order and is multiplied by the third 
power of the small parameter a/R. For typical experimental parameters mentioned above 
a/R ^ 10^'^. If one puts T = 1 K [the highest temperature where Eq. fl5Sl) is applicable] the 
T-dependent term on the right-hand side of Eq. fl38p appears to be of order 10"^^. 
Equation (138|) can be identically rewritten in the form 

^ / n^hcR ( a \ ^ ^ , ^. , . 

where the thermal correction to the Casimir force is given by 

It is interesting to compare this result obtained using the general formulation of the PFA 
with the exact result for the thermal correction in sphere-plate configuration jrj]. In the 
functional determinant representation the exact free energy of the Casimir interaction in 
sphere-plane geometry can be written as 

^a, T) = ^ ^ Tr In [1 - M(a, 6)] , (41) 

i=— oo 

where the explicit form of the matrix M in ideal metal case is contained in Ref. ^]. In 
the limit of extremely low temperature, T ^ Tr, and small separation distances, a ^ i?. 



Ref. 



74| arrives at the following exact result for the temperature-dependent part of the free 



-exact,. ( W ^29 , 112 a 



energy 

A.^r>.r) = ^^^f + (42) 

From this equation, the exact thermal correction to the Casimir force in the limiting case 
under consideration is the following 

^tK^ ia,T) = = - 1125 (her ■ ^^^^ 

It can be seen from the comparison of Eqs. (1401) and (1431) that although the most general 
formulation of the PFA gives the correct dependence of the thermal correction on the fourth 
power of T, the numerical coefficient before it is underestimated by a factor of 4.48. At 
the same time, the application of the simplified formulation of the PFA ( !T3|) in this limiting 
case results in an incorrect dependence of the thermal correction on T (the third power of 

16 



T instead of the fourth). These results are not surprising. As explained in Sec. II, the PFA 
is applicable only under the condition a ^ i? and leads to reliable predictions which are of 
zeroth order in the small parameter a/R. Because of this, any contribution either to the 
free energy or to the force which relative magnitude is numerically of about a/R or smaller 
calculated using the PFA must be disregarded as physically meaningless. Keeping in mind 
that in the limiting case of extremely low temperature, T ^ Tr, the relative magnitude of 
the thermal correction is much less than a/R, the thermal effect in the Casimir force in this 
case should be considered as unobservable. 

A completely different type of situation occurs under the condition T ^ Tr which is well 
satisfied in all experiments on measuring the Casimir force performed at room temperature 
(we remind that for the sphere of R = 100 yum radius Tr = 11.4 K). Here, irrespective of 
whether the inequality T ^ Tq or T ^ Ta is valid, the exact expression fHT]) leads in the 
zeroth order of the small parameter a/R to the simplified formulation of the PFA for the 

n 

thermal correction to the Casimir force between a sphere and a plate |74| : 

ArF37=*(a, T) = 27ri?AT J-pp(a, T). (44) 

In the case T ^ the thermal correction to the Casimir free energy for two plane parallel 
plates, ArJ^pp(a,T), is contained in Eq. f p7|) : 

Ar,.A.T,^-m(lX(^-l)- (45) 



16™^ \TaJ V Ta 

Thus, under a condition that T ^ the simplified formulation of the PFA is applicable 
not only to the total quantities Fsp(a,T) and J-pp(a, T) [see Eq. f|T3l) ]. but separately to the 
zero-temperature and thermal contributions to these quantities as well. This invalidates the 
statement of Ref. GgJ that for a sphere above a plate at T > Tr the PFA is inapplicable to 
each contribution alone made on the basis of the worldline approach. 

V. APPLICABILITY OF THE PROXIMITY FORCE APPROXIMATION TO 
DESCRIBE THERMAL CORRECTION TO THE CASIMIR FORCE BETWEEN 
REAL METALS 



Now we consider the possibility to apply the simplified formulation of the PFA for the 
theoretical description of the thermal correction to the Casimir force in recent experiments. 

17 



As noted in the Introduction, the thermal effect in the Casimir force is fundamentally dif- 
ferent depending on what model of the dielectric permittivity of metal (Drude or plasma) is 
used. It is most straightforward to illustrate this difference in the low-temperature regime, 
T ^ Ta, for two plane parallel plates described by the simple plasma model [Eq. flTBl) with 
no contribution of core electrons, i.e., with gj = 0] or the Drude model (1151) . In the case 
when the simple plasma model is used the temperature-dependent part of the free energy is 
given by [t^ 



he Jc(3) 



2 / 2^ ^ 4 



TaJ 90 \Ta 



So 
a 



(46) 




where 6q = 2'n'e/up is the effective skin depth in the frequency range of infrared optics. If, 



however, the Drude model is used under the condition T <^ one arrives at 



16| 



ArJ'(;\a,T)=ArJ'(l\a,T) 



(47) 



+ 



1 -4— + 12 

a \ a 



As is seen in Eq. ( H71) . the additional term emerging in the case of the Drude model is 
positive and linear in temperature. Thus, it dominates at low temperatures leading to an 
anomalously large thermal correction. Note that the condition of low temperature with 
respect to is satisfied in all modern experiments performed at room temperature at 
separations below 1 fim. At the same time for all these experiments the condition of high 



temperature with respect to Tr, i.e., T ^ Tpt, is satisfied. Because of this, it can be said [74 1 
that modern experiments on measuring the Casimir force and its gradient in a sphere-plate 
configuration belong to the region of medium temperatures. 

The measurement data of all experiments performed in the configuration of a rather large 



-55 



58 



with the theory using the PFA. In static experiments (e.g., in Refs . f^^.. .y^^ .y^ .^^ 
62]) Eq. ( fT3i) was used for this purpose. In dynamic experiments iJsjTjTjt^^ form 
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19^ 


isl,!: 


41: 





of the PFA in Eq. (!22i) was employed. As was explained in Sees. II and III, in both cases 
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terms of order a/i?, as compared with unity, were disregarded. This raises the question 
if this neglection imposes some constraints on the possibihty to observe thermal effects in 
the Casimir force. Below we analyze this question using both theoretical approaches to the 
calculation of thermal Casimir force discussed in Sec. II. 

We start with the Drude model approach using the tabulated optical data [68] for lme{u) 
extrapolated to low frequencies by means of Eq. (fT5|) . The thermal effect can be characterized 
by the relative thermal correction to the Casimir force 



Fsf\ a,T) - F^^ \a,0) 



4^'^) (a, T) = 'Zm^. Z". (48) 

and the Casimir pressure 



4-'(.r). ^'°''°^,-^^;;''°-») . (49, 

We compute the relative correction psj) by Eqs. ([5]) and f|T3l) over the separation region from 
0.1 to 5 /im. The relative correction fHOl) is computed over the same separation region by 
using Eqs. ([3]) and ([5]). We remind that the measurement data for the quantity P{a,T) in 
dynamic experiments are obtained by using Eq. (122|) from the gradient of the Casimir force 
dFsp{a,T)/da which in its turn was recalculated from the measured frequency shift. Thus, 
in static experiments using sphere-plate configuration, the PEA in the form of Eq. (fT3l) is 
part of the theory, whereas in dynamic experiments the PEA in the form of Eq. (122|) allows to 
convert the experimental data for the force gradient into the data for the equivalent Casimir 
pressure. 

In Fig. 5 we present the computational results in percents for the relative thermal cor- 
rections s!^'^^ (the solid line 1) and 5^'^'' (the solid line 2) at T = 300 K as a function of 
separation in the region (a) from 0.1 to 5 /im and (b) from 0.1 to 1 /im. In the same figure, 
the short-dashed and long-dashed lines show the quantity a/R m percents taken with the 
same sign as the thermal corrections as a function of separation for R = 100 and 150 /im, 
respectively. These lines demonstrate the size of typical errors arising from the use of the 
PEA. As can be seen in Eig. 5, the relative thermal corrections for both the Casimir force 
and equivalent Casimir pressure are rather large. In the region from 1 to 4 /im they take 
negative values and achieve -35% and -47%, respectively. Within the separation region 
from 0.1 to 1 /im the magnitudes of thermal corrections to the force and to the pressure 
increase from 1.5% and 0.7% to 23% and 16%, respectively. In the same separation region 

19 



the error due to the use of the PFA varies from 0.1% to 1% for the sphere of 100 fim radius 
and from 0.07% to 0.7% for the sphere of 150 /im radius. Thus, the use of the PFA in the 
Drude model approach does not impose any constraints on the possibihty to measure large 
thermal effects predicted in this approach. 

The answer to the question whether or not the thermal corrections to the Casimir force 
in a sphere-plate configuration can be measured depends not only on the errors due to the 
application of the PFA, but also on the size of the total experimental error. Thus, in the 
static experiment 55|] performed by means of an atomic force microscope the total relative 



experimental error of the Casimir force Fgp determined at a 60% confidence level at separa- 
tions 100 and 200 nm is equal to 3.5% and 22.4%, respectively These are larger than 
the re^ 



Refs. 



ative thermal correction d!^'^^ in Fig. 5(b). Because of this, in the experiment of 



the predicted thermal correction cannot be either confirmed or exc 



completely different type of situation occurs in dynamic experiments of Refs. 
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uded. A 



23| per- 



formed by means of a micromachined oscillator. Thus, in the most precise experiment of 



Refs. m 



231] the total relative experimental error of the Casimir pressure P determined 



even at a higher, 95%, confidence level varies from 0.19% to 0.9% and to 9% when sepa- 
ration increases from 160 to 400 and to 750 nm, respectively. As can be seen in Fig. 5(b), 
the magnitude of the relative thermal correction 5^'^^ remains much larger than the total 
experimental error over the entire range of experimental separations. That is the reason 
why, when no evidence of the predicted thermal corrections was found, Refs. 22, l23|] ar- 
rived at the conclusion that the Drude model approach is experimentally excluded at a 95% 
confidence level (within a narrower se par ation region from 210 to 620 nm this approach was 



excluded at a 99.9% confidence level 



23|). 



It has been speculated in the literature that the experimental exclusion of the Drude 
model approach might be not warranted because of some theoretical and experimental un- 
certainties. Specifically, it was stressed that the Drude parameters Up and -y used for the 
extrapolation of the optical data to low frequencies may vary. Besides, Ref. TSj surmised 
that measurements of absolute separations in Ref. (which data were also used to exclude 
the Drude model approach) contain an unaccounted systematic error which could bring the 



data away from the theory. Here, we make a test o 



to the theory-experiment comparison in Refs. 



22 



both these opportunities with respect 



23|. In Fig. 6 the differences between 



the theoretical Casimir pressures computed using the Drude model approach, as explained 



20 



above but with account of surface roughness, and the mean experimental Casimir pressures 



measured in Refs. 



22 



23| are indicated as dots. The theoretical pressures for the lower 



set of dots are computed with the conventional Drude parameters of Sec. II {up = 9.0 eV, 
7 = 0.035 eV). The theoretical Casimir pressures for the upper set of dots are computed with 
the alternative Drude parameters [3], oOp = 6.82 eV, 7 = 0.0405 eV, where the value of the 
plasma frequency is most different from the conventional one. The differences p*^<^°'' — pexpt 
with p^^^°^ computed using all other alternative Drude parameters listed in Ref. [73] are 
sandwiched between the two sets of dots shown in Fig. 6. The solid lines in Fig. 6 indi- 
cate the borders of the confidence intervals [— S(a),H(a)] computed in Ref. at a 95% 
confidence level. By the construction of [— S(a), S(a)], 95% of dots must belong to these 
intervals. However, as is seen in Fig. 6, almost all dots in both sets lie outside the solid 
lines demonstrating that the theoretical approach using the Drude model is excluded by the 
data at a 95% confidence level over the entire measurement range (as was mentioned above, 
within a bit narrower reparation region the exclusion at a 99.9% confidence level holds). 

n 

It is pertinent to note also that the use of optical data of Ref. [77[ for the first absorption 
bands (which are slightly different from the data of Ref. [68[ used in our computations) does 
not change this conclusion. 

Now we consider the possibility to bring the experimental data in agreement with the 
Drude model approach by assuming that there is an unaccounted systematic error in the 
determination of absolute separations. Note that this systematic error (if any) is separation- 
independent. This is because in the setup used the differences between the values of sep- 
arations where the Casimir pressures were measured are fixed interferometrically to high 
precision (see Ref. |20| for details). To find absolute separations, one should know the 
initial absolute separation which is determined from the electrostatic calibrations. If one 
admits that electrostatic calibrations contain some uncertainty, the initial separation would 
be burdened with some unaccounted constant systematic error which is translated to all 
separations. In an attempt to place the dots within the limits of the confidence intervals, 
in Fig. 7(a) we plot the differences pt^<=°'' — pcxpt -^j^j^ ^^^q gg^^g Qf ^j^g Drude parameters 
mentioned above as a function of separation, but with all separation distances decreased by 
Aa = 1 nm. (Note that the increase of separations results in even larger deviations of dots 
from the confidence intervals than in Fig. 6.) This corresponds to the assumption that an 



unaccounted systematic error in the measurement of separation distances in Refs. 



22 



23| is 
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equal to 1 nm (remind that in actual fact the total error in the measurement of separations 



in Refs. 



22 



23| is equal to 0.6 nm). As can be seen in Fig. 7(a), the shift of separations 



for 1 nm does not furnish the desired result for the lower set of dots at separations above 
230 nm and is not helpful for the upper set of dots at any separation. The decrease of all 
separations for 3 nm [see Fig. 7(b)] expels the lower set from the confidence intervals at 
separations below 230 nm, but not yet includes these dots into the confidence interval at 
a > 310 nm. In Fig. 7(c) we illustrate the largest decrease of separations for 6 nm. Here, all 
dots related to the lower set are still outside the confidence intervals at a > 430 nm, but the 
dots of both sets are already outside of them at shortest separations. Note that similar to 
Fig. 6 the differences pt^<=°'' — pcxpt computed with all other alternative Drude parameters 
are sandwiched between the two sets of dots shown in each of Figs. 7(a), 7(b), and 7(c). 
The same is correct when the optical data (not just the Drude parameters) from different 



sources are used. As shown in Ref. 
the data of Ref. 68 



79|, the use of some alternative optical data instead of 



would decrease the magnitudes of theoretical Casimir pressures and, 
thus, only increase discrepances between the predicted and experimental Casimir pressures. 
This means that no systematic error in the measurement of absolute separations can bring 



the experimental data of Refs. 



22 



23| in agreement with the Drude model approach using 



any values of Drude parameters and any sets of the optical data. 

From the above it follows that the sensitivity of dynamic experiments by means of mi- 
cromachined oscillator is quite sufficient to registrate the thermal correction to the Casimir 
pressure, as predicted by the Drude model approach. What is more, the application of the 
PFA in the form of Eq. ( l22l) in this case is warranted because both the total quantities 
dFgp/da and P and contained in them thermal corrections are much larger than a/R for the 
experimental parameters. The lack of any observation effect for the thermal correction to 
the Casimir pressure, as predicted by the Drude model approach, means that this approach 
is experimentally inconsistent. 

We are coming now to the question is it possible to observe small thermal effect in the 
Casimir force, as predicted by the plasma model approach in sphere-plate configuration. In 
this case the relative thermal correction to the force, and to the pressure, 6f^''^\ are 

expressed once again by Eqs. (HHj) and ( 149|) where Fsp^ and P*^'^-' are replaced with Fsp^ 
and P^^^ computed using the dielectric permittivity of the generalized plasma-like model 
( IT6l) . In Fig. 8 the computational results (in percents) are presented for 5^'^^ (the solid 
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line 1) and for 5^'^^ (the solid line 2) at T = 300 K as a function of a in the region (a) 
from 0.1 to 5 /im and (b) from 0.1 to 1 fim. The short-dashed and long-dashed lines show 
the quantity a/R (in percents) for R = 100 and 150 /im, respectively. As can be seen in 
Fig. 8, the thermal corrections computed using the plasma model approach are positive 
and increase monotonously with the increase of separation. At a < 1 /im their magnitudes 
are much less than the magnitudes of the thermal corrections computed using the Drude 
model approach (see Fig. 5). As a result, 6!^'^-' lies below the short-dashed line related to 
the static experiments. The relative correction to the Casimir pressure, lies below the 

long-dashed line related to the dynamic experiments at a < 1.6 /im. This means that one 
cannot attribute physical meaning to such small values of the thermal corrections 5^'^^ and 
5/^'^^ computed using the PFA [Eqs. ( JT3l) and (l22l) ] which disregards contributions to the 
force and to the pressure of order a/R. 

The magnitudes of the thermal corrections 5^'^^ and 5^'^^ are markedly larger than a/R 
shown by the short-dashed (long-dashed) lines only at separations a > 0.6 /im (a > 2/im). 
However, either in the static experiment IsS] on measuring the Casimir force Fsp or in 
the most precise dynamic experiment 22l, 1231] on measuring the Casimir pressure P, over 
the entire separations regions a > 0.1 /zm and a > 0.16 /xm, the magnitudes of 5^'^'' and 
J^'^-* fall far short of the total relative experimental error. It is not surprising, then, that 
thermal effects were not observed in these experiments. One more possibility for future 
measurements is to use spheres of larger radius. This allows to make a/R smaller and, thus, 
make thermal corrections 5^ '^^ and 5^'^^ computed using the PFA physically meaningful over 
the entire range of separations a > 0.1 /im. In addition, in static experiments the Casimir 
force Fsp is larger for a sphere of larger R resulting in a smaller relative experimental error. 
For spheres of centimeter-size radii, however, inavoidable deviations from perfect spherical 



shape is a problem of great concern which prevents accurate electrostatic calibration |80 |. 
It is therefore unlikely that small thermal effect in the Casimir force, as predicted by the 
plasma model approach, will be observed in the configuration of a sphere above a plate. 
We conclude this section with a brief discussion of recent exact results for the thermal 



Casimir force between a sphere and a plate described by simple 



plasma and Drude mod- 



els without account for interband transitions of core electrons [251] and by the generalized 
plasma-like and Drude- like models 26|. The key question is whether or not these results 
support computations performed using the PFA. As was discussed in Sec. IV, for ideal met- 
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als the exact results in the zeroth order in a/R coincide with the PFA if T ^ T/j. This is 
the case for both the total Casimir force or Casimir pressure and separately for the ther- 
mal corrections. Unfortunately, the exact computations of Refs. 25|, |26|] were performed 



in regions far away from the values of experimental parameters and outside the region 
where the PFA is applicable. For example, in Ref. 25| the sphere radii were chosen to be 
R = 0.1, 0.2, 0.5, 1, 2, 5 /im and computations were performed at separations from 0.5 to 
10 yum. For a sphere radius R = 10 fim the computations of Ref. 25| were made from a = 1 
to 10 /im. Thus, in all cases considered it holds a/R> 0.1. Remind that for the experimen- 



tal parameters of Refs 
Refs. 



55 



22 



76| a/R varies from 0.00063 to 0.003 and for the parameters of 



23| from 0.0011 to 0.005. 



According to the results presented in Sees. II and III, the PFA in the form of Eqs . ( IT^ 
and f l22|) is applicable only at a/R <^ 1. Because of this, the speculation of Ref. 25|] that 
at small separations the Drude and plasma models lead to Casimir force values much closer 
than predicted by the PFA made on the basis of computations in the region of a/R > 0.1 
is unjustified. The computational results of Ref. [25| at a = 0.5 /im, T = 300 K clearly 
demonstrate that deviations between the exact theory and the PFA vanish with decreasing 
a/R. Thus, from Fig. 3 of Ref. 25| it can be seen that relative deviation of the quantity 



Fsp (a, T) / Fs^j^"* (a, T) computed exactly from that computed using the PFA decreases from 
9.2% to 2.5% when the sphere radius increases from R = 0.1 fim {a/R = 5) to R = 5 fim 
{a/R = 0.1). In such a manner exact computations confirm that the computational error 
in the ratio Fsp\a,T)/Fs^\a,T) arising from using the PFA for real metals (0.025 for 
a/R = 0.1) is smaller that a/R. For the experimental values of a/R mentioned above this 
error would be much less than the total experimental error. Because of this it can be stated 
with certainty that the comparison of the performed experiments with the exact theory 
instead of the PFA (when the exact will become available) will not lead to any changes with 
respect to already obtained conclusions. 

Both papers 25|, |26| claim that for large separations in the sphere-plane geometry the 
Drude model leads to a force a factor of 3/2 smaller than the plasma model (instead of a 
factor of 2 as it holds for two parallel metal plates). In such a general form this formulation 
is, however, somewhat misleading as it does not cover all relevant limiting cases. In fact it 
is correct only in the case when Ap = 27Tc/ujp ^ i? ^ a |25|, |26|, i.e., outside the apphcation 
region of the PFA. If, however, we consider large separations a ^ hc/{2kBT), i.e., T ^ 
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Ta, and large sphere R ^ a (which means that the PEA is apphcable), then the exactly 



calculated ratio F}^'/F}p' = 2. The same result is obtained usine; the PFA. As shown in 



sp / ^sp 

Ref. |26| by means of exact computations, there is one more limiting case, a ^ R, R < Xp 
and T ^ Ta, where Fsp^ / Fsp^ = 1. This result is quite natural because a small sphere above 
aplate can be modelled as an atom characterized by some effective dynamic polarizability 



8l| . Keeping in mind that at large separations only the static polarizability is relevant which 



is common for the plasma and Drude models, both the free energies and forces computed 
using these models coincide in the limit of large a. Thus, in sphere-plate geometry the ratio 
Fsp^ /Fs^^ can take different values, 3/2, 2, and 1, in the large distance limit depending on 
the parameters of the problem. In the experimental situations, however, where a ^ i? and 
R ^ Ap, in the limit of large distances, a ^ hc/{2kBT), it holds Fsp'^ / Fsp^ = 2 like for two 
parallel plates. 
In Ref. 



26| the exact computations of the free energies ^s^p^ and ^sp ^ were performed 



for a sphere of i? = 5 /xm radius at T = 300 K and 77 K within the separation region from 
0.265 to 95 /im. This corresponds to the values of a/R varying from 0.053 to 19. Once 
again, only near the smallest a/R considered the PFA becomes applicable, whereas for the 
comparison with recent experiments more than one order of magnitude smaller ratios a/R 
are required. At the shortest separation a = 0.265 /im {a/R = 5.3%) from Fig. 5 of Ref. 26 1 
one obtains the exact value for the ratio J^sp/J^sp'^ = 1-11 at T = 300 K. The same ratio 
computed using the PFA is equal to 1.14 leading to a relative 2.6% deviation between the 
PFA and exact results which is less than one half of a/R. This confirms that for real metals 
in sphere-plate geometry the PFA results for both the thermal Casimir force and free energy 
go to the respective exact results with decreasing a/R. 



VI. THE POSSIBILITY TO MEASURE THERMAL EFFECT IN THE 
GRADIENT OF THE CASIMIR PRESSURE 

The original Casimir configuration of two parallel plates is the only one where the com- 
parison of experiment with theory does not require either the PFA or more sophisticated 
exact computational methods which can be effectively used so far solely in some restricted 
ranges of parameters. The only experiment of this kind was performed [82j] in the modern 
stage of the Casimir force measurements. This is an experiment of dynamic type where one 
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of the plates was oscillating at the natural oscillator frequency. As a result, the shift of 
this frequency under the influence of the Casimir pressure was measured and recalculated 
into the gradient of the Casimir pressure. The experiment 8^ reported 15% measure of 
agreement between the data and theory using the model of ideal metals (see review [ll| for 
details). 

The gradient of the Casimir pressure 

, dP(a,T) , , 

is obtained from Eqs. ([3]) and (|5]) by differentiation with respect to separation 

The thermal correction to the gradient of the Casimir pressure can be characterized by the 
quantity 

4.,,,).f:(-ii_z(M). (52) 

We will add additional indices D, p and IM and calculate the quantities 6!^'^\ 6!^'^^ and 
^(/M,3) ^gpgj^^jjj^g ^]^g used model of dielectric properties of metal described in Sec. II. 

The computational results for the relative thermal correction fl52|) to the gradient of 
the Casimir pressure as a function of separation at T = 300 K are presented in Fig. 9. 
Here, the separation region from 0.5 to 5 fim is considered because at shorter separations 
it would be hard to experimentally keep the plates parallel. The dashed, solid, and dotted 
lines represent results computed by using the tabulated optical data extrapolated to low 
frequencies by the Drude model with conventional parameters, the generalized plasma-like 
model, and the model of ideal metals, respectively. As can be seen in Fig. 9, for ideal 
metal plates and plates described by the generalized plasma-like model the relative thermal 
correction to the gradient of the Casimir pressure monotonously increases with the increase 
of separation. At separations a < 3 /im, where the experiment is feasible, 5^'^^ < 3.6% and 
^(/Af,3) ^ So small magnitudes of 5^'^^ make it impossible to observe the thermal 

corrections due to the generalized plasma-like model in the measurements of the gradient 
of the Casimir pressure. By contrast, the relative thermal correction 6!^'^^ computed using 
the Drude model approach decreases monotonously with increase of separation from 0.5 to 
4.1 fim where it achieves the minimum value of -53.3%. At a = 3 fim it holds 4^'^^ = -43%. 
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Thus, the proposed experiment on measuring the gradient of the Casimir pressure [83| might 
provide one more confirmation for the exclusion of the Drude model approach, but is not 
well adapted for the measurement of the thermal correction, as predicted by the generalized 
plasma-like model. Note that according to the authors of Ref. [ssl their effort to measure 
the Casimir force is not yet successful due to the use of Al surfaces. The point is that Al 
is the subject of quick oxidation even in high vacuum. For this reason, in the first precise 
measurement of the Casimir force jssi, Al surfaces were coated with thin transparent layers 
of Au/Pd. The use of such layers, however, complicates the comparison between experiment 
and theory. As a result, starting from 2000 almost all experiments on the Casimir force 



between metallic test bodies used Au coated surfaces. Presently another setup with 
parallel Au surfaces, instead of Al, is under construction jsS]. 

It can be shown also analytically that measurements of the gradient of the Casimir 
pressure are unsuitable for the detection of small thermal corrections. For example, for ideal 
metal plates in the limit of low temperatures T ^ it holds 



I [tV 120 T . ^ ij, 

3 \Tr, J TT T„ 



(53) 



Taking into account that the second contribution on the right-hand side of Eq. (1531) does 
not depend on separation, the main contribution to the gradient of the Casimir pressure in 
the low-temperature limit takes the form 

?^;^!^.gEi^^,,,-^^r,r-^. (54) 

oa 60a^ ^ ^ 

As can be seen on the right-hand side of Eq. fl54l) . there are no terms in powers of T/Ta 
and the thermal correction is exponentially small. Equation (15^ reproduces the exactly 
calculated dP^^^'^^ (a, T) /da at T = 300 K with an error of less than 1% up to the separation 
distance a = 3 fim. The absence of power-type thermal corrections in the gradient of the 
Casimir pressure for ideal metals explains why the thermal effect in this case is so small. 

At low temperatures, the analytical representation for the gradient of the Casimir pressure 
can be obtained for real metals as well. Thus, for metals described by the generalized plasma- 
like model at separations above 0.5 fim considered here the influence of interband transitions 
can be neglected. Then the thermal correction in the gradient of the Casimir pressure is 
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obtained as the negative second derivative of Eq. (H6!) leading to 



dP^P\a,T) 




) 



2 



da 



60a^ 



15C(3) (5o /T 
7r3 a \Ta 



(55) 



where the exponentially small terms are omitted. Here, the leading thermal correction is of 
power-type, but it contains the first and third power of small parameters. Equation fl55]) 



range of separations below 1.3 /im. Similar low-temperature results for metals described by 
the Drude model can be obtained by calculating the negative second derivative of Eq. fH7|) 



This approximate expression reproduces the values of dP^^^ (a, T) / da with an error less than 
1% at separations below 1 /im. 

To conclude this section, for the configuration of two parallel plates in the dynamic 
regime the relative thermal correction to the gradient of the Casimir pressure computed 
using the plasma model approach achieves 20% only at a = 5 /im (see Fig. 9). The same 
relative size of the thermal correction computed using the plasma model but to the Casimir 
pressure is achieved at a shorter separation a = 3.7 /im [see the solid line 2 Fig. 8(a)]. 
Thus, the static experiment with two parallel plates of sufficiently large area might be 
preferable for the registration of small thermal corrections at relatively large separations. 
Although at the same separation the relative thermal correction to the Casimir force in a 
sphere-plate configuration is larger [43% in accordance with the solid line 1 in Fig. 8(a) at 
a = 3.7/im], small values of forces for spheres of relatively small radii make its measurement 
impossible. As was discussed above, the use of centimeter-size spheres meets difficulties due 
to unavoidable deviations from sphericity. In this respect the increase of plate area and, 
thus, the increase of the Casimir force, may make possible to decrease the experimental 
error at separations of a few micrometers to about 10%, i.e., below the characteristic size of 
small thermal effect predicted by the plasma model approach. 



provides a better than 1% approximation for the exactly computed dP^P\a,T)/da over the 




(56) 
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VII. CONCLUSIONS AND DISCUSSION 



In the foregoing we have analyzed the possibihty to measure thermal effects in the Casimir 
force between metals in the configurations of a sphere above a plate and two parallel plates. 
In connection with the sphere-plate configuration, a detailed study of the reliability of the 
PEA to describe the thermal Casimir force was performed. We have considered the two 
formulations of the PFA, the most general, Derjaguin, formulation [49] and the simplified 



formulation of Ref. 



24( 1 which was used for the comparison of many experiments on the 



measuring of the Casimir force with theory. We show in Sec. II that for a large sphere, 
i.e., under the condition a <^ R, the magnitude of the relative difference between thermal 
Casimir forces Fsp computed for real metals by using the two formulations of the PFA is 
less than a/R. This holds regardless of what approach to the description of charge carriers 
(Drude or plasma model) is employed. A similar result is obtained in Sec. Ill for the gradient 
of the thermal Casimir force in sphere-plate configuration, dFsp/da. Here again the relative 
difference computed using the two formulations is less than a/R. The obtained results are 



compared with the experimental data of Ref. [72| demonstrating that for large spheres of 
different radii at T = 300 K the relative differences of the measured |-Fsp|; dFgp/da and those 
computed using the simplified formulation of the PFA are also less than a/R. From this 
we can conclude that the simplified formulation of the PFA is a good approximation for the 
calculation of thermal Casimir forces in sphere-plate configuration for spheres of large radii. 

In Sec. IV special attention is paid to the case of an ideal metal sphere above an ideal 
metal plate. In this case a simple analytic expression for the difference between the two 
formulations of the PFA is obtained. The thermal correction to the Casimir force computed 
using the PFA is considered both in the limiting case of very low temperatures, T ^ Tr, and 
in the case T ^ Tr. The latter includes the region of parameters where all measurements 
of the Casimir force in sphere-plate configuration have been performed. The results for the 
thermal correction obtained using the PFA are compared with available exact results. It 

n 

is emphasized that in the experimental region T ^ Tr the exact results of Ref. [7^ in the 
zeroth order of a/R coincide with the simplified formulation of the PFA. This provides a 
theoretical basis for the application of the PFA in the measurements of the thermal Casimir 
force for the comparison between experiment and theory. 

The applicability of the PFA and the possibility to observe the thermal correction to 
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the Casimir force in the configuration of a sphere above a plate made of real metals are 
discussed in Sec. V. Here, the relative thermal corrections to the Casimir force and to the 
Casimir pressure are computed in the framework of both the Drude model and the plasma 
model approaches. The obtained results are compared with the error introduced from the 
use of the PEA and with typical experimental errors of force and pressure measurements. 
New evidence is provided that some unaccounted systematic errors in the experiment of 



Refs. 
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23| cannot bring the measurement results in agreement with the prediction of the 



Drude model approach. It is argued that the configuration of a sphere above a plate is well 
suited to exclude the large thermal correction derived within the Drude model approach but 
scarcely can be used for the experimental observation of small thermal effects, as predicted 
from the plasma model approach. The reason is that the respective thermal correction 5^'' is 
several times larger than the error of PFA, a/ only at sufficiently large separations, where 
the error of force or equivalent pressure measurements exceeds 5'>^\ The exact computational 
results for the thermal correction to the Casimir force between real metals described by the 
Drude and plasma models are compared with the PFA results. It is shown that the region 
of parameters where the exact results are available is outside the area of application of the 
PFA. We demonstrate that at relatively large separations the ratio F^^/Fs^'^ can go to 



different limiting values (2, 3/2 or 1) depending on the re 
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ationship between the parameters 



26| . The performed comparison 



a, R, Xp and T. This generalizes the results of Refs. 
enables us to conclude that the available exact results for real metals converge to the PFA 
results when a/R ^ 0. For the experimental values oi a/R varying from 0.00063 to 0.005 
the errors due to the use of the PFA instead of the exact methods are much less than the 
experimental error of force and pressure measurements in all experiments performed to the 
present day. 

In the configuration of two parallel plates, we have calculated the relative thermal cor- 
rection to the gradient of the Casimir pressure which is measured in the dynamic regime 
(Sec. VI). We show that a small thermal effect, as predicted in the plasma model approach, 
is additionally suppressed in the gradient of the Casimir pressure. Because of this, the orig- 
inal Casimir configuration of two parallel plates of sufficiently large area used in the static 
regime remains the most prospective for measurement of small thermal effect in the Casimir 
force. 
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FIG. 1: The quantity / characterizing differences between the two formulations of the PFA for the 
Casimir force as a function of (a) separation at room temperature T = 300 K and (b) temperature 
at a separation a = 2 fim. The dotted, dashed and sohd Unes show the results computed using 
ideal metals, Drude and plasma model approaches, respectively. The sphere radius R = 100 /xm. 
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FIG. 2: The quantity / characterizing differences between the two formulations of the PFA for 
the Casimir force as a function of a/R at room temperature T = 300 K. The dotted, dashed and 
sohd lines show the results computed using ideal metals, Drudc and plasma model approaches, 
respectively. The groups of lines numerated 1, 2, and 3 are for the respective fixed separations 
a = 0.1, 2, and 5 /xm. 
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FIG. 3: The quantity p characterizing differences between the two formulations of the PFA for the 
gradient of the Casimir force as a function of (a) separation at room temperature T = 300 K and 
(b) temperature at a separation a = 2 jim.. The dotted, dashed and sohd hues show the results 
computed using ideal metals, Drude and plasma model approaches, respectively. The sphere radius 
R = 100 /xm. 
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FIG. 4: The quantity p characterizing differences between the two formulations of the PFA for 
the gradient of the Casimir force as a function of a/R at room temperature T = 300 K. The 
dotted, daslicd and solid lines show the results computed using ideal metals, Drude and plasma 
model approaches, respectively. The groups of lines numerated 1 and 2 are for the respective fixed 
separations a = 0.1 and 5 /xm. 
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FIG. 5: Relative thermal correction to the Casimir force (the solid line 1) and to the Casimir 
pressure (the solid line 2) computed using the Drude model approach at T = 300 K over the 
separation region (a) from 0.1 to 5/xm and (b) from 0.1 to 1/Ltm. The short-dashed and long- 
dashed lines show the quantity a/R for R = 100 and 150 /xm, respectively. 
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FIG. 6: Differences between tfie tfieoretical Casimir pressures, computed by means of the Drude 
model approach, and mean experimental Casimir pressures at different separations are shown as 
dots. The lower and upper sets of dots correspond to the use of conventional and alternative Drude 
parameters (see text for further details). The solid lines indicate the borders of the confidence 
intervals found at a 95% confidence level. 
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FIG. 7: Differences between the theoretical Casimir pressures, computed by means of the Drude 
model approach, and mean experimental Casimir pressures at different separations are shown as 
dots. Compared to Fig. 6, all separation distances are decreased by (a) Aa = 1 nm, (b) Ao = 3nm 
and (c) Aa = 6 nm. 
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FIG. 8: Relative thermal correction to the Casimir force (the solid line 1) and to the Casimir 
pressure (the solid line 2) computed using the plasma model approach at T = 300 K over the 
separation region (a) from 0.1 to 5/um and (b) from 0.1 to 1 /xm. The short-dashed and long- 
dashed lines show the quantity a/R for R = 100 and 150 //m, respectively. 
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FIG. 9: Relative thermal correction to the gradient of the Casimir pressure between two parallel 
plates as a function of separation computed at T = 300 K using ideal metals (the dotted line) , and 
the Drude and plasma model approaches (the dashed and solid lines, respectively). 
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